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Abstract 

A  set  of  DSMC  computations  using  the 
SOCRATES  computer  code  has  been  performed 
to  aid  in  the  development  and  analysis  of  the  MirEx 
experiment  to  observe  UV  radiation  from  space¬ 
craft  thruster  plumes.  The  plumes  of  interest  for 
this  report  originate  from  the  amine  propellant  main 
engine  of  the  Progress-M  and  Soyuz-TM  space¬ 
craft.  The  altitude  of  the  plume  observations  is 
approximately  380  km.  The  SOCRATES  calcula¬ 
tions  provide  a  baseline  estimate  of  the  spatial  dis¬ 
tribution  of  the  exhaust  gases  as  they  interact  with 
the  ambient  low-density  atmosphere  composed 
predominantly  of  atomic  oxygen.  In  addition,  a  pro¬ 
posed  mechanism  for  OH(A)  excitation  based 
upon  high-velocity  collisions  between  H20  and  O  is 
exercised.  A  sequence  of  calculations  was  per¬ 
formed  to  investigate  the  influence  of  the  angle  of 
attack  of  the  engine  upon  the  predicted  radiation 
levels.  In  general,  it  is  found  that  the  predicted 
OH(A)  generation  has  a  very  large  spatial  scale,  on 
the  order  of  several  kilometers.  Comparison  of  the 
predictions  with  data  showed  the  predictions  were 
high,  leading  to  speculation  that  the  excitation 
cross  section  used  in  the  calculations  may  be  sus¬ 
pect.  A  preliminary  variation  of  the  energy  thresh¬ 
old  was  performed  using  a  set  of  experimental  con¬ 
ditions.  It  was  found  that  an  increase  in  the  energy 
threshold  could  bring  the  predictions  into  agree¬ 
ment  with  both  the  spatial  distribution  and  the 
absolute  signal  level  of  the  experimental  data. 

Introduction 

Exoatmospheric  rocket  exhaust  plumes  present 
a  challenging  situation  for  accurate  modeling 
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because  the  rarefied  reactive  atmosphere  interacts 
with  the  plume  over  very  large  length  scales. 
Understanding  the  dynamics  of  exhaust  plumes  or 
jets  is  critical  to  missions  such  as  docking  maneu¬ 
vers  near  future  space  stations,  where  the  plume 
may  impinge  on  critical  surfaces.  In  addition  to 
understanding  the  gas  dynamics,  the  post-com¬ 
bustion  reactivity  of  plumes  is  also  of  interest.  In 
particular,  it  has  been  reported  that  visible  and 
ultraviolet  glows  exist  in  the  plume  far  downstream 
of  the  nozzle  exit.  Recently,  an  experimental  effort 
was  initiated  to  measure  the  ultraviolet  radiation 
from  very  high  altitude  rocket  exhaust  plumes.1,2 
The  concept  is  to  use  the  Mir  space  station  as  a 
measurement  platform  to  conduct  passive  radio- 
metric  observations  of  the  rocket  exhaust  jets  of 
spacecraft  which  supply  the  Mir  station.  As  stated 
in  Refs.  1  and  2,  an  ultraviolet  imager  aboard  Mir 
has  recorded  radiation  emitted  from  plumes  in  the 
bandpass  from  280  to  320  nm.  Because  this  radia¬ 
tion  is  robust  and  easily  detected,  the  experimental 
effort  has  attempted  to  identify  the  radiator  and  iso¬ 
late  the  chemical  mechanism  responsible  for  the 
excitation.  This  is  intended  to  be  the  first  step  of  a 
multitask  process  that  will  eventually  include  data 
acquisition  on  many  of  the  UV  emitters  in  other 
bandpasses. 

As  discussed  in  Ref.  2,  it  is  surmised  that  the 
detected  radiation  originates  from  the  OH(A) 
specie,  although  some  contribution  to  the  radiation 
may  be  due  to  NH(A).  Viereck,  et  al.,  have  reported 
considerable  NH(A-X)  radiation  in  the  plume  of  the 
small  VRCS  US  Shuttle  motors,  while  both  OH(A- 
X)  and  NH(A-X)  have  been  identified  in  the  larger 
PRCS  thruster.3  Since  the  instrument  response 
curve  of  the  current  Mir  UV  instrumentation  heavily 


*  The  research  reported  herein  was  performed  by  the  Arnold  Engineering  Development  Center  (AEDC),  Air  Force  Materiel 
Command.  Work  and  analysis  for  this  research  were  performed  by  personnel  of  Sverdrup  Technology,  Inc.,  AEDC  Group,  technical 
services  contractor  for  AEDC.  Further  reproduction  is  authorized  to  satisfy  needs  of  the  U.  S.  Government. 

t  Senior  Member,  AIAA. 

*  Member,  aiaa.  Approved  for  public  release;  distribution  unlimited. 

This  paper  is  declared  a  work  of  the  U.  S.  government  and 

not  subject  to  copyright  protection  in  the  United  States.  ^ 

American  Institute  of  Aeronautics  and  Astronautics 


AIAA-99-0975 


emphasizes  OH(A-X)  radiation  at  310  nm  over 
NH(A-X)  radiation  at  337  nm,  it  will  be  presumed 
for  the  purposes  of  this  numerical  study  that  the 
dominant  radiation  source  measured  in  the  MirEx 
imagery  is  OH(A-X). 

The  authors  have  previously  put  forth  a  hypoth¬ 
esis  for  the  excitation  of  OH(A),  originally  due  to 
Kofsky,  et  al.4  The  hypothesis  is  that  H20  mole¬ 
cules  are  expelled  from  the  orbiting  nozzle  at  a  rel¬ 
atively  high  velocity  with  respect  to  the  “quiescent” 
atmosphere.  Since  orbital  velocities  are  on  the 
order  of  7  to  8  km/sec,  and  typical  plume  nozzle 
exit  velocities  can  be  on  the  order  of  3  km/sec,  a 
considerable  relative  velocity  may  exist  between 
the  plume  molecules  and  the  atmospheric  atomic 
oxygen.  The  largest  relative  velocity  of  the  plume 
H20  with  the  atmosphere  is  the  case  in  which  the 
motor  is  retrofiring.  Kofsky,  et  al.,  reviewed  the  lab¬ 
oratory  measurements  for  H20  +  O  ->  2  OH(X)  and 
reasoned  that  the  production  rate  of  OH(A)  would 
be  proportional  to  the  ground-state  OH  production 
rate  by  a  factor  depending  upon  an  activation 
energy.  They  estimated  the  activation  energy  as 
simply  the  heat  of  formation  required  to  produce 
OH  in  the  excited  electronic  state.  In  terms  of 
experimental  parameters,  the  threshold  velocity  to 
overcome  the  assigned  activation  energy  is  10.45 
km/sec.  While  this  is  certainly  a  reasonable 
approach,  no  experimental  data  have  been 
recorded  on  the  exact  form  of  the  reaction  cross 
section  for  direct  OH(A)  production  from  high- 
energy  collisions  between  H20  and  O. 

The  purpose  of  the  present  Direct  Simulation 
Monte  Carlo  (DSMC)  modeling  of  rocket  exhaust 
plumes  is  twofold  (1)  to  provide  guidance  for 
experiment  planning,  and  (2)  to  assess  the  likeli¬ 
hood  of  possible  excitation  schemes.  The  resolu¬ 
tion  of  the  problem  will  require  an  iterative  method¬ 
ology  between  experiment  and  theory  until  the 
point  is  reached  when  the  excitation  and  the  reac¬ 
tion  scheme  can  be  unequivocally  identified  and 
simulated.  Certainly,  the  combined  experimental 
and  theoretical  effort  has  not  reached  this  point 
yet,  as  there  are  many  issues  that  are  yet  to  be 
resolved.  In  this  context,  it  is  still  very  fruitful  to 
learn  the  characteristics  of  the  H20  +  O  reaction 
mechanism,  and  to  assess  the  consistency  of  the 
data  acquired  to  date  with  the  numerical  results. 


Clearly,  if  the  numerical  results  cannot  be  recon¬ 
ciled  with  the  data,  then  a  different  hypothesis  for 
OH(A)  excitation  should  be  sought. 

For  this  analysis,  two  data  collection  activities 
will  be  evaluated.  Both  the  Progress-M  and  Soyuz- 
TM  spacecraft  used  the  same  rocket  motor, 
denoted  in  this  report  as  the  PME.  In  one  case,  the 
Progress-M  underwent  dedicated  maneuvers,  and 
data  were  acquired  at  near  90-deg  aspect  angle  to 
the  plume  flow.  In  the  second  case,  the  Soyuz-TM 
crew  transport  undertook  standard  reentry  proce¬ 
dures  for  retrofiring,  accomplishing  the  4-min  dura¬ 
tion  principal  burn  in  darkness.  While  the  analysis 
will  demonstrate  properties  of  the  interaction  of  H20 
and  the  ambient  atomic  oxygen,  the  purpose  is  to 
understand  the  spatial  lengths,  energy  scales,  and 
collision  dynamics  of  the  flow  in  general.  As  an  aid 
to  future  experiment  planning,  two  numerical  varia¬ 
tions  are  performed,  one  in  which  the  angle  of 
attack  is  varied,  and  the  other  as  a  variation  of  the 
molecular  energy  barrier  for  OH(A)  production. 
These  variations  are  intended  to  assess  the  H20  + 
O  hypothesis  and  determine  the  sensitivity  to  exper¬ 
imental  conditions.  As  a  last  action,  the  most 
recently  acquired  data  will  be  used  to  compare  with 
the  absolute  level  of  predicted  OH(A)  generation  as 
a  test  of  reasonableness  of  the  H20  +  O  hypothesis. 

Method 

As  mentioned  above,  the  DSMC  results  pre¬ 
sented  here  will  focus  on  replicating,  to  an  extent, 
two  measurements  of  the  PME  exhaust  plume. 
While  much  greater  detail  may  be  found  in  Refs.  1 
and  2,  a  brief  description  is  given  here.  The  PME 
is  the  main  thruster  used  in  the  Progress-M  cargo 
ship  and  the  Soyuz-TM  crew  transport  ship.  The 
engine  uses  UDMH  and  N204  as  fuel  and  oxidizer, 
respectively,  with  an  O/F  ratio  of  1 .8.  The  nominal 
thrust  of  the  PME  is  3100  N  (700  Ibf).  It  has  a  noz¬ 
zle  expansion  area  ratio  of  48.  Measurements  of 
the  near-field  exhaust  plume  radiation  indicate  the 
PME  has  a  rapid  startup.2  The  PME  is  quite  often 
used  to  boost  the  Mir  orbit  while  the  Progress-M  is 
docked.  A  numerical  solution  for  the  combustion 
and  nozzle  flow  has  been  reported  earlier,1  and 
the  nominal  exit  plane  characteristics  are  shown  in 
Table  1. 
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Table  1.  Computed  Exit  Plane  Gas  Flow 
Parameters  of  the  PME  Rocket 
Engine 


Temperature,  K 

698 

Pressure,  N/m2 

1.05x103 

Density,  molecules/cm5 

1.09  xIO17 

Velocity,  cm/sec 

2.99  x  10s 

[H20] 

0.2884 

[N2] 

0.2668 

[H2] 

0.1905 

[CO] 

0.1891 

[C02] 

0.0531 

[H] 

0.0110 

[OH] 

i  .07  x  i  o-5 

[O] 

9.20  x  ICT6 

]  mole  fractions 


In  the  first  measurement,  which  occurred  on  15 
May  98,  an  unmanned  Progress-M  spacecraft  was 
positioned  near  the  Mir  station  for  the  execution  of 
three,  short  duration  (5-15  sec)  burns.  Two  of  the 
burns,  denoted  as  dV3  and  dV4,  were  performed 
such  that  the  Progress-M  plume  was  within  the 
field  of  view  of  the  UV  imager  located  inside  the  Mir 
station.  The  details  of  the  burns  are  summarized  in 
Table  2. 


Table  2.  Summary  of  PME  Retrofirings 


Reentry 

Spacecraft 

Progress-M 

Date 

1 5  May  98 

15  May  98 

Duration,  sec 

15 

5 

244 

Range  from  Mir  to  PME,  km 

00 

9,2 

15  to  28 

Imager  viewing  angle,  deg 

103 

84 

173  to  168 

Thrust-Wind  angle,  deg 

157 

176 

180  to  164 

PME  Altitude,  km 

386 

381 

385  to  377 

The  second  key  measurement  was  performed 
on  August  25,  1998,  while  the  Soyuz-TM  manned 
spacecraft  began  its  initial  descent  to  earth  after 
having  been  docked  to  Mir  on  a  10-day  mission. 
The  Soyuz-TM  craft  dedocked  from  Mir  and  drifted 
to  a  position  15  km  behind  Mir  and  at  2-km  higher 
altitude.  The  PME  was  ignited  at  that  point.  The 
retrofiring  was  lengthy,  lasting  over  240  sec.  As  the 
burn  occurred,  the  Soyuz-TM  slowed  slightly  and 
decreased  altitude  from  385  km  to  377  km.  In  addi¬ 
tion,  the  plume  thrust  vector  gradually  deviated 
from  the  vehicle  velocity  vector  due  to  orbital 


mechanics  at  a  rate  of  roughly  0.066  deg  per  sec¬ 
ond.  Thus,  the  firing  began  with  a  180-deg  angle 
between  the  plume  thrust  vector  and  the  wind  vec¬ 
tor,  and  finished  with  a  164-deg  angle  between  the 
plume  thrust  vector  and  the  wind  vector. 

The  DSMC  code  used  to  model  the  PME  flow 
was  SOCRATES,  written  by  Elgin,  et  al.,  in  1992 
for  the  USAF  Research  Laboratory,  formerly 
known  as  Phillips  Laboratory.5  The  SOCRATES 
code  was  originally  written  to  study  plume  contami¬ 
nation  effects  on  the  U.  S.  Shuttle  produced  by 
burns  of  the  Shuttle  motors.  The  code  has  been 
successfully  used  in  the  analysis  of  the  NH(A-X) 
emission  reported  by  Viereck  in  Ref.  3. 

The  underlying  method  of  SOCRATES  is  based 
upon  the  work  of  Bird.6  The  adaptation  provided  by 
Elgin,  et  al.,  incorporates  a  nonuniform  three- 
dimensional  Cartesian  grid.  The  code  input  allows 
the  user  to  specify  the  desired  chemistry  set  of 
interest  completely.  The  chemical  processes  can 
be  input  either  in  Arrhenius  form  or  as  tabulations 
of  energy  and  cross  section.  For  Arrhenius  rates, 
SOCRATES  performs  a  standard  manipulation  to 
determine  the  cross  section  as  a  function  of  the  rel¬ 
ative  velocity  of  the  collision  partners,  using  the 
variable  hard  sphere  model.  The  internal  energy  of 
the  molecules  is  taken  into  account  by  the  Larsen- 
Borgnakke  phenomenological  model.7 

Rather  than  perform  a  detailed  continuum  fluid 
dynamical  analysis  of  the  nozzle  exit  flow,  the 
SOCRATES  code  replaces  the  nozzle  with  a  point 
source  and  uses  Brook’s  model  to  determine  the 
distribution  of  molecules  entering  the  solution 
domain.8  This  is  entirely  justified,  since  the  solution 
domain  has  length  scales  of  kilometers,  while  the 
nozzle  diameter  is  only  0.3  m. 

In  the  current  application,  a  computational 
domain  was  defined  which  encompassed  a  volume 
of  10  km  x  10  km  x  13  km  about  the  PME  in  the  x, 
y,  and  z  dimensions,  respectively.  The  plume 
thrust  axis  is  placed  along  the  +z  dimension,  with 
the  wind  vector  placed  along  the  -z  dimension  for 
a  full  retrofiring.  The  domain  was  divided  into  7168 
cells,  using  a  scheme  of  16  x  16  x  28  grid  planes. 
The  SOCRATES  default  value  of  three  simulated 
molecules  per  minor  species  per  cell,  and  20  sirnu- 
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lated  molecules  per  major  species  per  cell  was 
used,  although  a  variational  study  of  these  values 
did  not  appreciably  affect  the  results.  (A  major  spe¬ 
cies  is  defined  to  have  a  mole  fraction  within  a  fac¬ 
tor  of  three  of  the  maximum  mole  fraction  in  a 
given  cell.) 

SOCRATES  determined  the  atmospheric 
parameters  used  in  the  current  calculation,  via  an 
implementation  of  the  MSIS  '90  atmospheric 
model.  For  example,  at  an  altitude  of  383  km,  the 
free-stream  number  density  was  determined  to  be 
1 .34  x  1 08  molecules-cm-3.  The  ambient  tempera¬ 
ture  was  883  K,  which  results  in  a  mean  free  path 
of  the  undisturbed  atmosphere  as  2.4  x  106  cm, 
and  a  nominal  collision  rate  for  a  given  molecule  of 
1  collision  per  22  secs.  The  molar  composition  of 
the  atmosphere  at  this  altitude  is  approximately  94- 
percent  atomic  oxygen,  4-percent  molecular  nitro¬ 
gen,  and  2-percent  helium,  although  these  values 
depend  upon  the  solar  cycle,  etc. 

For  the  present  modeling,  the  exhaust  flow  was 
assumed  to  be  nonreactive,  with  the  exception  of 
the  proposed  H20  +  O  reaction.  The  lack  of  a 
detailed  reactive  model  for  the  near-field  plume 
exhaust  gases  prohibits  a  close  examination  of  the 
near-field  flow,  on  the  order  of  tens  of  meters,  as 
far  as  OH(A)  generation  is  concerned.  That  exami¬ 
nation  is  left  to  a  later  work.  For  now,  the  primary 
task  is  to  estimate  the  far-field  structure  of  the 
plume  and  assess  the  plausibility  of  the  H20  +  O 
mechanism  for  OH(A)  generation. 

Results 

As  a  first  step,  the  PME  flow  for  a  nominal  retro- 
firing  was  calculated.  The  PME  exit  plane  proper¬ 
ties  (Table  1)  were  used  at  an  altitude  of  383  km. 
The  PME  was  assumed  to  have  a  vehicle  velocity 
of  7.35  km/sec,  in  keeping  with  the  experimental 
value.  The  plume  was  directed  directly  into  the 
oncoming  atmospheric  wind,  which  is  composed 
mostly  of  atomic  oxygen. 

The  predicted  total  number  density  is  shown  in 
Fig.  1  as  a  cut  through  the  solution  domain  along 
the  symmetry  axis.  The  dimensions  of  the  solution 
region  are  13  km  horizontally  and  10  km  vertically. 
The  PME  nozzle  is  located  5  km  from  the  left  side 


of  the  solution  region,  and  is  denoted  by  the  '+.' 
The  grid  size  near  the  nozzle  exit  is  240  m  in  the 
axial  direction  and  380  m  in  the  other  two  direc¬ 
tions.  Thus,  the  rapid  expansion  of  the  plume  gas 
from  the  initial  density  of  1017  cm-3  is  not  fully  cap¬ 
tured  by  this  coarse  resolution.  The  ambient  den¬ 
sity  is  1.3  x  108  cm-3.  The  outer  contour  in  Fig.  1, 
marked  A,  indicates  a  density  of  2  x  108  cm-3, 
while  the  contour  line  marked  B  indicates  a  density 
of  2.5  x  109  cm-3.  The  B  contour  is  roughly  2  km 
downstream  of  the  nozzle  exit  and  provides  a  rea¬ 
sonable  length  scale  for  the  far-field  plume. 


13  km 


Fig.  1.  Total  number  density  in  the  PME  retrofiring 
plume. 

Figure  2  shows  a  profile  of  the  total  number 
density  taken  along  the  plume  axis  of  Fig.  1 .  Again 
we  note  that  the  ambient  density  is  ~108  cm-3;  the 
lower  values  of  the  number  density  "behind"  the 
PME  indicate  that  the  PME  plume  "sweeps"  the 
ambient  atmosphere,  creating  a  spatial  region  of 
lower  density.  Also  shown  is  the  gradual  decay  of 
the  far-field  plume  density  to  the  ambient  atmo¬ 
spheric  density. 


Fig.  2.  Axial  profile  of  the  total  number  density  in 
the  PME  retrofiring  plume. 
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The  density  of  H20  exhausted  in  the  plume  will 
be  proportional  to  the  total  density,  since  H20  is 
modeled  as  comprising  a  fixed  proportion  of  almost 
30  percent  of  the  plume.  Fig.  3  displays  the  mean 
velocity  of  the  expelled  water  in  each  computa¬ 
tional  cell.  The  position  of  the  PME  is  denoted  in 
the  figure  by  a  "+."  The  mean  velocity  is  obtained 
by  averaging  over  all  H20  molecules  in  a  given 
cell,  including  the  molecules  that  have  undergone 
scattering  collisions,  as  well  as  those  that  have  not. 
In  this  context,  the  mean  velocity  provides  a  picture 
of  the  dominant  behavior  of  the  molecules.  It  must 
be  emphasized,  though,  that  the  average  velocity 
of  a  given  species  can  be  a  misleading  parameter 
when  speaking  of  single-collision  events  such  as 
the  production  of  OH(A)  from  the  high-energy  colli¬ 
sion  of  an  H20  molecule  with  an  oxygen  atom.  The 
direction  and  magnitude  of  the  mean  velocity  vec¬ 
tor  should  not  be  interpreted  as  precluding  individ¬ 
ual  molecular  trajectories  and  collision  events.  The 
data  for  Fig.  3  were  extracted  from  a  slice  through 
the  center  of  the  solution  domain.  The  length  of  the 
vectors  indicates  the  relative  scaling  of  the  mean 
velocity,  while  the  color  of  the  vectors  corresponds 
to  the  z  component  of  the  mean  velocity  vector. 
The  plume  exhaust  is,  in  general,  turned  back  into 
the  direction  of  the  atmospheric  wind,  but  over  dis¬ 
tance  scales  of  tens  of  kilometers.  The  z  velocity 
component  along  the  plume  axis  is  plotted  in  Fig. 
4.  The  vacuum  acceleration  of  the  nozzle  exhaust 
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Fig.  3.  H20  average  flow  velocity  in  PME  retrofiring 
plume,  colored  by  average  Z-component. 


flow  is  noted  by  the  increase  in  the  H20  velocity 
from  the  nozzle  exit  value  of  3  x  1 05  cm/sec  to  over 
3.2  xIO5  cm/sec  roughly  1  km  downstream. 


Fig.  4.  Axial  profile  of  the  H20  average  Z-compo- 
nent  of  the  flow  velocity  in  the  PME 
retrofiring  plume. 


An  important  parameter  is  the  ambient  atomic 
oxygen  distribution,  taken  in  the  reference  frame  of 
the  spacecraft.  Figure  5  shows  the  density  con¬ 
tours  of  the  O-atom  for  the  center  slice  through  the 
solution  domain,  where  the  PME  is  denoted  by  the 
"+."  The  most  notable  feature  is  the  depletion  of  O- 
atoms  behind  the  spacecraft,  where  even  after  5 
km  from  the  PME,  the  O-atom  hole  shows  no  sign 
of  filling.  In  front  of  the  PME  one  can  discern  a 
slight  rise  in  the  O-atom  density,  presumably  due 
to  a  "snowplow"  effect.  It  was  found  that  the  plume 
exhaust  does  not  significantly  alter  the  mean  atmo¬ 
spheric  flow  direction.  Of  course,  the  individual 
direction  of  travel  of  those  O-atoms  which  scatter 


Fig.  5.  O-atom  number  density  in  the  PME  retrofir¬ 
ing  plume. 
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off  the  plume  molecules  can  be  radically  altered  as 
a  result  of  the  collision. 

The  collision  frequency  of  atomic  oxygen  and 
H20  is  shown  in  Fig.  6.  It  is  somewhat  deceptive  in 
that  the  contour  map  resembles  the  total  density 
contours  of  Fig  1.  However,  the  location  of  the 
PME  is  marked  on  Fig.  6  by  the  "+."  This  figure 
shows  that  the  maximum  in  the  0-H20  collision 
frequency  occurs  nearly  1  km  downstream  of  the 
PME  nozzle  exit. 


13  km 


Fig.  6.  Collision  frequency  of  H20  with  O-atoms  in 
the  PME  retrofiring  plume. 


In  addition  to  the  180-deg  thrust-wind  angle,  a 
90-deg  thrust-wind  angle  was  to  estimate  the  spa¬ 
tial  scale  of  the  flow  when  the  exhaust  is  being 
blown  sideways.  Figure  7  presents  the  computed 
H20  density  contours,  which  will  be  nearly  identical 
to  the  total  density  contours.  Figure  7  shows  that 
the  rarefied  atmosphere  does  not  significantly  dis- 


Fig.  7.  H20  number  density  in  the  PME  perpendic¬ 
ular-firing  plume. 


turb  the  flow  direction  of  the  PME  exhaust,  at  least 
to  within  the  first  5  km  available  in  the  computation. 
However,  the  H20  mean  flow  velocity  is  notably 
distorted  from  the  180-deg  solution,  as  shown  in 
Fig.  8.  In  this  case,  the  leeward  side  of  the  PME 
exhaust  expands  quite  uniformly  but  the  windward 
side  of  the  mean  flow  velocity  is  not  a  uniform 
expansion.  The  color-coding  of  Fig.  8  indicates  the 
z-velocity  component  for  each  computational  cell 
and,  as  expected,  only  the  portion  of  the  exhaust 
most  windward  to  the  atmosphere  will  have  suffi¬ 
cient  relative  velocity,  on  average,  to  overcome  the 
OH(A)  reaction  threshold. 


13km 


Fig.  8.  H20  average  flow  velocity  in  PME  perpen¬ 
dicular-firing  plume,  colored  by  average  Z- 
component. 


The  atomic  oxygen  density  is  shown  in  Fig.  9 
for  the  90-deg  thrust-wind  angle  case.  In  this  con¬ 
tour  plot,  the  deficit  of  O  atoms  behind  the  PME 
plume  is  again  apparent  as  in  the  180-deg  case  of 
Fig.  5,  but  with  the  variation  that  the  deficit  region 
is  stretched  to  shadow  the  entire  plume.  Also 
apparent  in  Fig.  9  is  the  excess  density  of  O  atoms 
on  the  windward  side  of  the  plume.  The  peak  of  the 
excess  is  roughly  a  factor  of  two  over  the  ambient 
level.  An  interesting  point  is  that  the  O-atom 
excess  is  not  located  at  the  nozzle  exit,  where  the 
plume  density  is  highest.  Rather,  it  is  located 
roughly  1.2  km  from  the  PME  nozzle  exit,  presum¬ 
ably  where  the  optimal  combination  of  density  and 
spatial  divergence  creates  the  most  effective  "net" 
for  collecting  the  fast  O  atoms.  This  is  displayed  in 
the  map  of  the  O-atom  mean  velocity  shown  in  Fig. 
10.  The  color  coding  again  refers  to  the  z-compo- 
nent  of  the  mean  velocity,  and  it  is  seen  that  in  the 
region  downwind  of  the  plume,  the  O-atoms  are 
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significantly  slowed,  and  very  slightly  "entrained" 
by  the  exhaust  flow,  resulting  in  an  upward  curva¬ 
ture  to  the  flow  direction.  The  direction  straightens 
very  quickly,  however,  although  the  mean  O-atom 
velocity  downstream  of  the  plume  is  much  slower 
than  the  ambient  value. 


13  km 


Fig.  9.  O-atom  number  density  in  the  PME  per¬ 
pendicular-firing  plume. 


13km 


Fig.  10.  O-atom  average  flow  velocity  in  PME  per¬ 
pendicular-firing  plume,  colored  by 
average  Z-component. 


The  examination  of  the  radiation  from  the  PME 
exhaust  flow  is  a  closer  comparison  to  the  measur¬ 
able  data  than  the  gasdynamic  calculations. 
Although  SOCRATES  does  not  contain  a  radiative 
model  for  OH(A),  it  does  come  with  a  utility  to  con¬ 


struct  line-of-sight  summations  of  the  flow  proper¬ 
ties.  This  tool  is  usually  quite  useful  for  predicting 
various  emitters  in  the  IR,  such  as  C02,  for  which 
radiative  models  exist  within  SOCRATES.  In  this 
case,  however,  the  volumetric  generation  of  OH(A) 
will  be  sufficient  to  replicate  the  expected  spatial 
distribution  of  OH(A)  radiation.  Since  the  lifetime  of 
OH  (A)  is  very  short,  6.9  x  10-7  sec,  the  transport 
and  quenching  of  OH(A)  can  be  neglected,  and  the 
rate  of  reactions  which  produce  OH(A)  can  be 
treated  as  equivalent  to  the  rate  at  which  OH(A) 
photons  are  emitted.  Under  these  assumptions,  it 
is  possible  to  construct  a  two-dimensional  spatial 
maps  of  OH(A)  emission  in  equivalent  units  of  pho- 
tons-sec-1-cm-2. 

Figure  11  is  an  OH(A)  emission  map,  in  units  of 
photops-sec_1-cm-2  obtained  by  integrating  the 
solution  domain  along  parallel  lines  of  sight.  The 
OH(A)  excitation  mechanism  was  assumed  to  be 
the  H20  +  O  reaction,  with  an  activation  energy  of 
110  kcal/mol  (equivalent  to  10.45  x  105  cm/sec). 
This  computation  encompassed  a  solution  space 
of  10  km  x  10  km  x  20  km,  and  the  PME  nozzle  exit 
is  shown  in  the  figure  for  clarity.  The  general  shape 
of  the  emission  matches  the  expectations,  given 
the  divergence  of  the  plume  density  shown  in  Fig. 
1 .  Note  that  the  peak  intensity  occurs  on  the  plume 
centerline,  roughly  2  km  downstream  of  the  nozzle. 
This  is  quite  easily  seen  on  the  axial  profile  of  Fig. 
12.  In  that  plot,  the  axial  distance  is  referenced  to 
the  PME  nozzle  exit.  The  OH(A)  intensity  rises 
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Fig.  11.  OH(A)  emission  map  produced  with  base¬ 
line  assumptions. 
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over  the  first  2  km  of  plume.  Downstream  of  2  km, 
the  plume  intensity  gradually  declines  in  the  aver¬ 
age,  but  with  notable  oscillations.  Further  review 
with  finer  grid  spacing  will  determine  if  the  oscilla¬ 
tions  are  an  artifact,  or  a  true  predicted  structure. 
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Axial  Distance,  m  x  100 

Fig.  12.  Axial  profile  of  the  projected  OH(A)  emis¬ 
sion  map. 

A  key  experiment  variable  is  the  emission  inten¬ 
sity  as  a  function  of  angle  between  the  thrust  vec¬ 
tor  and  the  wind  vector.  Figure  13  is  a  composite  of 
ten  different  SOCRATES  computations  in  which 
the  angle  between  the  thrust  vector  and  the  wind 
vector  decreased  from  180  to  90  degs  in  10-deg 
increments.  For  each  computation,  the  OH(A) 
emission  map  was  generated  from  the  line-of-sight 
projections  of  the  solution  domain.  In  Fig.  13,  each 
OH(A)  map  has  dimensions  of  14  km  x  14  km, 
although  the  solution  domain  only  occupies  13  km 
x  10  km  of  that  region.  Furthermore,  the  plume 
thrust  angle  is  rotating  in  successive  images,  such 
that  the  wind  vector  for  all  the  images  is  horizon¬ 


tally  from  the  right  to  left.  The  angles  on  the  top  row 
decrease  from  180  deg  to  140  deg,  while  the  bot¬ 
tom  row  images  decrease  from  130  deg  to  90  deg. 

In  the  first  two  increments,  the  general  shape 
and  scale  of  the  OH(A)  emission  is  relatively 
unchanged,  even  though  the  plume  thrust  angle 
has  moved  from  180  deg  to  170  and  160  deg.  The 
intensity  of  the  emission  declines  as  the  thrust 
angle  is  rotated  to  150  degrees,  and  further  to  140 
degrees.  However,  notice  that  the  line  of  maximum 
emission  is  still  close  to  the  wind  vector,  and  not  the 
plume  thrust  vector.  This  is  in  keeping  with  the  fairly 
high  relative  velocity  requirement  for  OH(A)  pro¬ 
duction  in  the  H20  +  O  mechanism.  As  the  thrust 
axis  is  rotated  even  further,  we  see  a  significant 
decline  in  the  OH(A)  emission,  attributed  to  the 
decline  in  the  off-axis  number  density  of  H20 
exhausting  from  the  nozzle  (see  Fig.  1),  and  com¬ 
pounded  by  the  angular  variation  of  the  exit  flow 
velocity  from  the  maximum  along  the  thrust  axis 
(see  Fig.  3).  Finally,  for  the  thrust  vector  perpendic¬ 
ular  to  the  wind  vector,  there  is  only  a  minor  amount 
of  OH(A)  production,  generated  solely  on  the  wind¬ 
ward  side  of  the  plume.  The  total  plume  intensity  is 
shown  as  a  function  of  the  thrust-wind  vector  angle 
in  Fig.  14.  The  decline  in  the  total  OH(A)  intensity, 
as  measured  in  the  number  of  reactions-sec'1  gen¬ 
erating  OH(A),  declines  by  over  a  factor  of  30  as 
the  thrust  vector  rotates  from  retrofiring  to  perpen¬ 
dicular  with  respect  to  the  wind  vector  due  to  the 
change  in  available  translational  energy  to  over¬ 
come  the  activation  energy  requirement. 


a.  180  deg  b.  170  deg  c.  160  deg  d.  150  deg  e.  140  deg 
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Fig.  13.  Effect  of  angular  orientation  on  OH(A)  emission  for  varying  angles  af  attack. 
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Fig.  14.  Variation  in  total  OH(A)  produced  as  a 
function  of  the  angle  between  the  thrust 
and  wind  vector. 


The  preceding  computations  were  performed 
using  the  values  that  Kofsky,  et  al.,  had  set  forth 
concerning  the  activation  energy  required  for  the 
OH(A)  generation  to  proceed.  The  reaction  rate 
used  to  this  point  has  been, 

H20  +  O  ->  OH(X)  +  OH(A) 

koH(A)  =  3.8  x  1CT15  T1-3  e('Ea/kT) 

with  Ea  =  110.4  kcal/mol  (4.79  eV).  The  prefactor 
value  for  the  rate  and  the  temperature-dependent 
factor  are  based  upon  a  reporting  of  ground-state 
OH  production  by  Cohen  and  Westberg.10 

In  an  attempt  to  investigate  the  sensitivity  of  the 
production  of  OH(A)  to  this  energy  threshold,  com¬ 
putations  were  done  in  which  the  threshold  energy 
was  varied,  while  the  angle  of  attack  was  fixed  at 
180  deg.  Fig.  15  shows  a  composite  of  the  results 
for  the  OH(A)  emission  maps.  As  the  activation 
energy  increases,  in  this  composite,  the  location  of 
the  intensity  maximum  moves  closer  to  the  nozzle 
exit  as  the  overall  intensity  declines.  Note  that 
each  (a),  (b),  and  (c)  have  their  own  relative  color 
map  scales,  in  units  of  reactions-sec"1  -cm'2.  The 
relative  change  in  the  activation  energy  for  each 
adjacent  image  is  0.1  eV  in  (a),  0.2  eV  in  (b),  and 
0.4  eV  in  (c).  A  decrease  in  absolute  intensity  of 
the  OH(A)  with  increasing  Ea  is  observed. 

Figure  16  summarizes  the  activation  energy 
computations  by  performing  the  spatial  integrals  of 
the  OH(A)  emission  maps  of  Fig.  15.  In  this  figure, 
the  sharp  decline  in  the  emission  of  OH(A)  for 
increasing  Ea  is  clearly  shown.  It  is  quite  interesting 
in  Fig.  16  to  observe  the  small  rise  in  intensity  at 
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c.  0.4  eV  increments 

Fig.  15.  Variation  in  the  OH(A)  emission  as  a  func¬ 
tion  of  activation  energy. 


Fig.  16.  OH(A)  production  as  a  function  of  variation 
in  the  H20  +  O  reaction  activation  energy. 


roughly  Ea  =  5.1  eV.  Note  that  several  additional 
computations  were  done  about  this  feature  in  order 
to  resolve  it.  One  explanation  for  this  peak  is  that  as 
the  activation  energy  increases,  the  spatial  location 
of  the  intensity  maximum  slowly  moves  toward  the 
nozzle.  At  5.1eV,  an  optimum  combination  of  H20 
density  and  O-atom  density  combines  to  give  an 
increase  in  the  predicted  OH(A)  intensity.  Of 
course,  the  activation  energy  is  a  numerical  not  a 
physical,  variable,  and  the  appearance  of  the  peak 
is  not  an  experimentally  useful  feature.  Nonethe¬ 
less,  the  OH(A)  emission  level  is  very  sensitive  to 
the  value  chosen  for  the  activation  energy.  For 
example,  a  10-percent  upward  adjustment  of  the 
activation  energy  from  4.79  eV  will  result  in  a  60- 
percent  reduction  in  the  OH(A)  produced  in  the  flow. 
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At  this  point,  a  preliminary  comparison  to  the 
data  of  August  25,  1998  can  be  done.  Using  the 
baseline  computation  conditions,  i.e.,  Ea  =  4.79  eV, 
and  the  thrust-wind  vector  angle  =  174  deg,  a  pre¬ 
dicted  OH(A)  distribution  can  be  obtained  which 
corresponds  to  the  geometry  of  the  measurement. 
The  OH  (A)  emission  map  is  shown  in  Fig  17,  in 
which  the  intensity  is  in  photons-sec-1  -cm-2.  The 
computational  result  of  Fig.  17  has  spatial  dimen¬ 
sions  of  10  km  x  10  km.  The  image  data  acquired 
from  the  Mir  station  on  August  25, 1998  at  the  start 
of  the  Soyuz-TM  reentry  burn  have  an  angular 
dimension  of  8  deg,  at  a  range  from  the  Mir  of  15 
km.  A  typical  data  image  is  shown  in  Fig.  18  for 
comparison  with  the  prediction  of  Fig.  17.  The  data 
image  has  been  enhanced  so  that  the  bright  region 
corresponding  to  the  intense  near-field  radiation  is 
saturated,  while  the  far-field  luminous  region  is  evi¬ 
dent  across  the  circular,  active  portion  of  the  image. 
Note  that  a  rough  estimate  of  the  spatial  dimension 
of  the  active  portion  of  the  image,  i.e.  the  circular 
region  within  the  image,  has  a  diameter  on  the  order 
of  2  km.  As  noted  in  much  more  detail  in  Ref.  2,  the 
estimated  power  emitted  by  the  PME  plume  within 
the  UV  imager  field  of  view  is  -180  W.  Note  that  the 
contribution  of  the  near-field  radiation  to  the  total 
value  is  less  than  a  percent.  The  radiated  power 
corresponds  to  an  emission  of  approximately  3 
x  1020  photons/sec.  Restricting  the  spatial  integra¬ 
tion  of  the  predicted  plume  shown  in  Fig.  17  to  an 
area  roughly  2  km  in  diameter  results  in  a  predicted 
value  of  3  x  1 022  photons-sec-1 ,  which  is  two  orders 
of  magnitude  higher  than  the  observation.  Even 
allowing  for  a  reasonably  high  uncertainty  in  the 
measured  data,  the  current  reaction  mechanism  will 
overestimate  the  radiation  either  due  to  a  cross  sec¬ 
tion  which  is  too  large,  or  failures  in  the  model  such 
as  collisional  quenching  of  OH(A),  or  its  implemen¬ 
tation,  such  as  inappropriate  grid  cell  sizes,  to  rep¬ 
licate  the  true  flow  character. 

One  speculative  explanation  for  the  overpredic¬ 
tion  is  the  presence  of  an  energy  barrier  to  the 
reaction.  As  noted  above,  the  production  of  OH(A) 
is  rather  sensitive  to  the  value  used.  For  example, 
the  solution  that  had  an  activation  energy  of  6.8  eV 
was  manipulated  to  produce  a  map  with  the  August 
25.  1998  measurement  viewing  parameters.  The 
result  was  that  the  predicted  emission  was  2  x  1 020 
OH(A)  emitted  photons-sec-1,  which  is  in  fair 
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Fig.  17.  Predicted  OH(A)  emission  for  August  25, 
1998  Soyuz-TM  reentry  burn,  with  Ea  =  4.8 
eV. 


2  km 


Fig.  18.  Image  of  Soyuz-TM  PME  plume  acquired 
from  Mir  Station  during  August  25,  1998 
PME  operation. 

agreement  with  the  data,  considering  the  rough 
order  of  magnitude  approximations  used. 

Discussion 

This  analysis  has  attempted  to  provide  the 
framework  for  experimental  data  analysis  and 
future  experiment  planning.  The  global  flow  param¬ 
eters  were  examined,  characterized  by  H2O  for  the 
plume  exhaust.  The  rarefied  atmosphere,  com¬ 
posed  almost  totally  of  atomic  oxygen  (as  pre¬ 
dicted  by  the  MSIS'90  model),  permits  the  PME 
exhaust  plume  to  extend  over  a  distance  several 
kilometers  long.  The  PME  disturbs  the  ambient 
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environment  to  a  great  extent,  causing  the  deple¬ 
tion  of  the  atomic  oxygen  behind  the  plume  by  an 
order  of  magnitude.  This  was  true  regardless  of  the 
orientation  of  the  plume  with  respect  to  the  wind 
vector,  as  demonstrated  in  calculations  of  ram  and 
perpendicular  firings. 

A  predicted  OH(A)  emission  map  was  gener¬ 
ated  for  a  retrofiring  PME  on  the  basis  of  the 
hypothesized  H2O  +  O  reaction  mechanism.  That 
result  indicated  the  maximal  plume  radiation 
occurs  downstream  of  the  nozzle  exit  by  1  to  2  km. 
This  is  consistent  with  the  previous  preliminary 
modeling  done  at  350  km  altitude  in  Ref.  1.  How¬ 
ever,  the  predicted  far-field  radiation  is  higher  than 
the  measured  data  by  a  considerble  amount. 

The  overprediction  of  the  OH(A)  emission  using 
the  baseline  activation  energy  is  not  inconsistent 
with  the  H20  +  O  mechanism.  An  underprediction 
of  the  radiation  would  cast  doubt  upon  the  mecha¬ 
nism,  since  the  SOCRATES  model  is  an  idealistic 
case  as  it  lacks  potential  quenching  and  radiative 
effects.  While  these  effects  are  expected  to  be 
small,  they  will  lessen  the  prediction.  Ideally,  a 
closer  result  would  be  desirable;  however,  the  over¬ 
prediction  is  likely  due  to  an  overestimate  of  the 
reaction  cross  section,  i.e.,  the  rate  constant.  The 
H20  +  O  reaction  equation  contains  three  parame¬ 
ters  in  the  expression  of  the  rate  constant;  the  pref¬ 
actor,  the  temperature  exponent,  and  the  activation 
energy.  As  a  first  look  into  the  effect  of  uncertainty 
in  the  rate  constant,  an  assessment  of  the  OH(A) 
production  with  varying  activation  energy  was  done 
for  a  retrofiring  PME  plume.  Recall  the  origin  of  the 
rate  constant  expression  as  being  the  electronic 
ground-state  reaction  rate  expression,  with  the  acti¬ 
vation  energy  replaced  by  the  enthalpy  of  the 
excited  state  products.  Assuming  that  the  prefactor 
and  the  temperature  exponent  describe  the  mea¬ 
sured  dynamics  of  the  H20  +  O  ground-state  colli¬ 
sion  system,  it  is  not  unreasonable  to  consider  the 
activation  energy  as  a  source  of  uncertainty  for  the 
electronically  excited  system. 

Computations  of  the  electronically  excited  inter¬ 
action  hypersurfaces  of  H20  and  O  may  be  used  to 
confirm  or  reject  the  concept  of  an  activation 
energy  barrier  to  the  production  of  OH(A)  and  OH. 
Furthermore,  such  a  computation  would  be  able  to 


provide  a  prediction  on  the  other  factors  in  the  rate 
constant  expression,  as  well  as  insight  on  the 
expected  form  and  distribution  of  the  internal 
energy  of  OH(A).  This  could  then  be  tested  against 
calibrated  spectra  obtained  in  future  MirEx  tests. 
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